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Near  Axis  Ship  Wakes 


NEAR  AXIS  SHIP  WAKES 


I.  Introduction 


The  problem  we  address  Is  that  of  the  structure  of  a  ship 

surface  wave  wake  far  down  stream  in  the  vicinity  of  the  axis, 
t  y 

(i.e.:  -  ,  -  <  <  1,  where  L  is  a  typical  dimension  of  the  ship). 

Basic  approximations  made  are: 

(i)  The  linearized  theory  is  used. 

2 

(ii)  It  is  assumed  that  g  L/U  >  >  1.  Further,  rather 
conventional  stationary  phase  methods  are  employed  It  Is 
indicated  that  this  is  the  most  suspect  part  of  the 
calculation.  Elucidation  is  warranted. 

We  are  particularly  interested  In  the  effects  of  finite  ship 
width.  Accordingly,  we  deal  with  a  simplified  hull  model.  It  is  one 
of  a  prolate  semi-spheroid.  This  has  the  advantage  that  the 
underlying  protentlal  flow  is  simple  and  can  be  written  in  closed 


The  main  conclusions  obtained  are: 


1)  In  agreement  with  other  calculations  it  is  found  that 
singularities  appear  in  the  flow  field.  However,  the 
singularity  is  not  on  the  axis  but  is  displaced  therefrom. 
(In  our  model  the  displacement  is  of  the  order  of  the  ship 
width.) 

2)  It  is  argued  that  the  singularity  is  not  real,  but  is 
rather  due  to  an  injudicious  rise  of  the  stationary  phase 
method.  However,  it  is  probable  that  the  singularity 
indicates  that  large  disturbances  can  result  far  from  the 
Kelvin  angle  and  yet  not  on  the  axis. 
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II.  The  Basic  Equations 


I.  The  Integral  Representation 

These  are  as  In  reference  (1).  We  use  the  linearized  Euler 
equations  and  boundary  conditions.  In  a  coordinate  system  at  rest 
with  respect  to  the  ship  we  can  write  the  velocity  potential  as 

*  -  Ux  +  4>  (1) 

(U  is  the  velocity  of  the  ship  which  is  moving  in  the  +  x 
direction). 


Then  $  satisfies 


V2 


♦ 


0. 


(2) 


and 


3n 


-n  U  (on  the  ship  surface) 


(3) 


— £  +  ■  0  (on  the  free  ocean  surface  z  •  0)  (4) 

3x2  U2 


If  we  introduce  a  Green's  function  satisfying 


7  G<r'.  r)  -  6(r'  -  r) 


2  2 
U  3 


—  <r.  £>  +i~  -^r  G  <£'*  £> 

37'  »  Vt**" 


Z'-  0 


then  using  Green's  identity  we  have  the  representation: 


MO  -  /  G(r',  t)  -gp-  0  (£',  r)  dSj 


+  r  c-  /  *d')  ijr  «  <£'•  O  -  #  0  <!'.  £) 


Z'-O 

X'-X-(y') 


g  c+  J  ^  ax'  G  (£  *  V  ax'  G  *  £) 


arG(£'’  £>  dy' 

z  -o 

X'»X+(y') 


(Here  S_  is  the  hull  of  the  ship  and  c+  +  c-  is  the 


waterline.) 


The  equation  (7)  plays  two  roles. 


a)  When  r  is  on  this  is  an  inhomogeneous  integral  for 


4>(r)  . 


b)  When  r  is  not  on  this  gives  an  integral 
representation  for  $(r)  in  terms  of  its  value  and  normal  derivative 
on  S  j . 
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Our  basic  approximation  Is  then  the  following: 


where 


with 


and 


where 

expect  G 
r  on 

equation 

problem- 


The  G  defined  by  equation  (5)  and  (6)  can  be  written  as 

G  E  )  “  G0  (£'»  E>  +  Gi  <£'•  E>  (8) 


G 

o 


D 


{ 

i  r-i 


+ 


lr-r| 


(9) 


£  *  (*,  y.  -z). 


•  E) 


2  2 

i  d  k  k  ,  .  , 

-1  rr  _ X  k(z+z  ) 

<2.)2  V) 

u 


[k  (x'-x)+k  (y'-y)] 
e  x  y 

(10) 


2 

When  gL/U  >  >  1  and  r,  r'  are  on  ,  we 
<  <  Gq  .  Thus  the  Integral  equation  (7)  for  <p(r)  with 
can  be  approximated  by  putting  G  ■  G  .  The  Integral 
which  results  is  then  that  of  a  simple  potential  flow 
•which  is  readily  solved  yielding  a  solution  $ 


To  calculate  the  radiated  field  (for  r  far  from  S^)  we  then 
replace  $  In  equation  (7)  by  and  G  by  G^. 

More  precisely,  we  use  the  radiative  part  of  G^.  This  is 
obtained  so:  Causality  tells  us  that  the  poles  in  equation  (10)  are 
to  be  above  the  contour  of  integration.  The  radiative  part  is  just 
the  contribution  of  the  poles.  For  the  part  of  G^  which  generates  the 
radiation  field  we  then  obtain 


i 

2w 


-8.  fa 
2  J-a 

u 


d*y  Zx  y2  Z(z+z') 
1-2  Z  ® 


[£x(x'-x)+iy(y'-y)] 


(ID 


+  complex  conjugate 


(Here  we  have  introduced  dimensionless  variables 

so  that  (k  ,  k  ,  k)  -  -f  (Z  ,  Z  ,  l )  .  )  (12) 

x  y  yi  x  y 

2.  Boat  Model 


We  model  the  hull  as  half  of  a  prolate  spheroid.  To  be 


specific,  introduce  spheroidal  coordinates  so  that 


y  -  c  sinh  n  sin  9  cos  u> 


Z  -  c  sinh  n  sin  0  sin  a) 


Then  the  hull  is  the  surface 


n  **  n  ,  0  <  0  <  it,  tt  <  to  <  2tt 

o  —  —  —  — 


(The  major  and  minor  axis  are  thus  a  *  c  cosh  b  ■  c  sinh  nQ 


The  scale  factors  are 


-  2  2  2 

hj  ■  h2  *  c  /  sinh  n  cos  0  +  cosh  n  sin  0 


h^  =»  c  sinh  n  sin  0, 


or  if  we  introduce  £  -  cosh  n,  v  *  cos  0 


.  .  /  2  2 
hl“h2*c/?  *  U 


h,  «  c  /  (C2-  1)  (l  -  m2)  • 
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(2 


dS,  =*  c  /  c  —  u  /  c  -  1  sin  0  d0  du  . 
i  o  o 


Let  us  see  how  some  typical  integrals  are  to  be  evaluated 


Consider  1^.  This  is  an  Integral  over  c+.  c+ 
is  X'  >  0,  Z'  -  0 


0  <  0  <  - 


For  y'  >  0  ,  we  have  y'  ■  b  sin  0 


dy"  »  b  cos  0  d  0 


As  y  goes  from  0  to  b,  0  goes  from  0  to  jr  . 

2 


For  y'  <  0,  we  have  y'  *  -b  sin  0,  dy'  ■  -  b  cos  0  d  0 


As  y'  goes  from  -b  to  0,  9  goes  from  t  to  0. 


Clearly  in  our  approximation  of  ~  >  >  1  ,  the  principal 

U2 

contributions  come  from  the  vicinity  of  u>  *  x,  2x  . 


Thus, 


l  b  sin  0  »x 

p  y 


J  *  e 


U 


— £  Jib  sin  9  sin  w 

.2  .  U 

/  du>  e 


-Ml  Jt  b  sin  8  ,2x  .  — §  Jib  sin  0  sin 


+  e  y2  y 


3V  dto>  e  u2 


In  the  first  of  these  integrals,  we  let 


(i)  -  ir  +  e,  then  do>  -  de 


sin  w  ■  sin  (x  +  e)  -  cos  x  sin  e  *  -e. 

u2 

gJl  b  sin  8 

Similarly, 

2 

-2x  w  U 

jx^  gJl  b  sin  9 

2 

and  thus 


Then 


■  l 


de  c 


-gi  b  sin  9  e 
U2 
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u 


gi  b  sin  8 


-lib  sin  8  l  b  sin  8 

r  1  U  .  V  l 

t  e  +  e  }  • 


Inserting  Into  equations  (27)  and  (26)  yields: 


I,  -  <n>  If  <77>  / 


-  it  t  ^  1  + 

x_  X  U 


3  *"  3?  v2ir'  *  t  (1-2 1) 

“CO 

i  — ^  [£x  a  cos  8  +  b  sin  8] 


/n  cos  0  d8  {  e  U 


(29) 


i  — *•  (l  a  cos  8  -l  b  sin  8] 

,  ir  x  y  x 

+  e  1 


The  remaining  1^  are  evaluated  similarly.  The  results  are: 


,  .  a  “i“o  (1  x  +  l  y] 

....  d  l  l  ,,2  x  y 

I  +l  +1  .  /“  --X.-5  e  U 

12  13  2*  1-21 


rtr  u (C  -1) 
/  coa  0  d0  (  -^  [ 
o 


§i 

3C  ‘  2  2  iJ  ’  2  2 


1  <0*Il-»] 

-  Tl  + 


«o 


( V  "u  > 


{  e++  e"  }  .  (30) 


+  c.  c. 


and 


At  i 


-i-f  (t  X  +  tv) 


IV.  The  Stationary  Phase  Evaluation 

So  far  the  approximations  made  are: 

(i)  Linearization 

(ii)  $  +  ^n  S^.  This  justified  by  the  assumption 
gL/U2  >  >  1  . 

(iii)  In  the  integrals  over  only  the  immediate  vicinity 
of  the  surface  contributes  significantly.  Again  this  is  justified  by 

gL/U2  >  >  1  . 

We  are  left  with  generic  integrals  of  the  form 

-i-f  (i  x  +  l  y) 

«  TT  * 

I  +  -  /  dl  g(i  )  e 

—  — oo  y  ° '  y' 

i-f-  fi,  a  cos  8  ±  &  b  sin  0  1 
X  h(0)  e  u  d  0  (32) 

(Remember:  i  -  /”T  ,  4  -  /"i2  +  i2  .  ) 

x  x  y 

We  want  to  evaluate  these  integrals  in  the  limits 
ly/xl  <  <  1,  |^/x|  <  <  1  •  Our  approach  is  to  use  the  method  of 


stationary  phase  twice. 


First  we  do  the  9  Integral  by  stationary  phase.  The 
stationary  points  6-  are  then  functions  of  ly.  The  ly  Integral  la 
then  done  by  another  application  of  stationary  phase. 


The  main  result  is  that  in  the  limit  of  interest  there  are 
two  classes  of  stationary  points. 


(a)  4  ~  4  -  1,  i  ~  0, 

x  y 


cos  9  *1 

o 


(b)  4  ~  4  >  >  4  >  >  1,  cos  0-0 

y  x  o 


In  the  Integrals  I  —  ,  we  first  encounter 


I  h(6)  e1  *£ 


where  [4  a  cos  0  +  4  b  sin  0} 


.2  1  x 


3*+ 


The  stationary  phase  condition  ^  -  0 


At  • 

yields  tan  9 ,  m  +  — 


(33) 


1  $+  (9+) 


h(9+)  / 


/  h(0)  e1  ■  e 


-  ~T  (8-frf)2^  ( 9+) 
ez  d0 


Then 


-  -  4  (6-9  )2  *+(0+)  -i* 

I +  -  J  g(*y)  h(0+)  t  /  e  -  d 9 ]  e  —  (34) 


where 


*  «  — [1  (x  -  a  cos  9+)  +  1  (y  +  b  sin  0  +] 

*  IT  x  y 


(35) 


34 

The  stationary  phase  condition  —  ■  0 

o  X 


yields 


+  b  sin  0  +  -  y  +  1  b  cos  0  +  -  1  a  sin  0  +  30+ 

x-  .  -  y  -x  -  . 


3*x 

31  x  -  a  cos  0  +  +  ^ 

y 


x  -  a  cos  0  + 


]  ^  •  06) 


We  are  interested  only  in  cases  where  |x|  >  >  a.  Therefore, 
we  are  always  justified  in  taking  as  our  basic  equation 


31  +b  sin  6  +  -  y  +1  b  cos  0+  -  l  a  sin  0  +  30+ 

X  —  —  J  —  V  —  X  —  — 


31 


+  [- 


31 


(37) 


Since  0+  are  implicit  functions  of  1^  via  equation  (33) 
this  is  in  principle  a  very  difficult  equation  to  solve  for  1^  . 
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However,  in  our  limit  it  is  readily  solved  by  a  perturbation  theory. 
We  assume  the  second  term  on  the  right  is  small  compared  to  the 
first.  Calculate  4^  using  only  the  first  term  and  then  verify  that 
the  second  term  is  indeed  small. 


Thus,  in  first  approximation  our  equation  is 


34 

> 

34_ 


(38) 


where 


±b  sin  0  -  y 

o 


,  and  by  assumption 


I  0  +  I  <  <  I  • 


Since 


3* 


34  4  (24-1) 

y  * 


Our  basic  equation  (38)  becomes 


4  (24-1) 
x 


(39) 


Squaring  gives 


4X2  (24-1 ) 


(40) 
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£  and  set  21-1  -A,  we 


If  we  note  that  £  2  «  £(£-1),  £  2* 

y  x 


obtain  a  quadratic  equation  for  A  which  can  be  solved  to  give 

1  +  /l  -  8  ei 


(41) 


4  e; 


Clearly,  the  two  solutions  correspond  to  the  wave  trains 
whose  caustics  give  the  classical  Kelvin  19  cone. 

First,  we  look  at  the  -  in  equation  (41).  This  gives  us  the 
simpliest  (and  least  interesting)  of  the  two  wave  sets.  For  small 
0  this  gives 

A  -  1-2  Q2  . 


Then 


£  -  -  !  -  02 


(42) 


and 


£  .  /I  „  i  ° 

x  2 


(43) 


From  equation  (39)  we  then  see  that 


£  -  Q  (44) 

v  +  ' 
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that 


What  does  this  imply  for  0q  ?  The  equation  (33)  tells  us 


tan  0  +  *  +  0  —  (45) 

—  —  a 

Thus  9+  is  close  to  either  0  or  it  depending  on  the  sign  of 
the  right  side  of  equation  (45). 

Then  0  +  ■  -y/x  .  (Here  x  <  0  and  we  will  always  look  at  y  >  0  .) 
Then  9+  -  0  »  -y/x 

9_  -  tt  -  e 

and  solving  for  e  we  obtain 

£  =  0  *  -y/x  ,  i.e.  9  »_x  -  0 

Then  *+  (0+)  -  -  . 

From  this  it  follows  that 

i  <e-e+)2  ♦+«*>  +  i  ,/4 

d9  -  /  -  e- 

ga 


f  e 


The  equation  for  the  6+  are 


tan  6+  =  ±  I  w. 


To  be  definite,  let  us  discuss  the  case 


y  >  b.  Then  0+  >  0  . 


0+  (0_)  is  then  just  less  than  (greater  than) 


Explicitly, 


0 


+ 


20+a 

_ 


20  a 


Substituting  in  the  expressions  for  <{>+(&+)  yields 


w 


h/2 


and  then 


00 


4  ( 6-0+)  <K  (9+) 


I  e 
—00 

y/®rrU2  e£  + 


d0 


iir/4 


gb 


1+ 


/ °°^U2  0+  g+i^A  e-i$+  (  0+,  0+) 
gb  -  ~  - 


-f  (A  -A  ) 
2  y  y' 


2 

2  3  4 


3A  dA 

y  y 


Here 


4+  -  -f  [Ax  -  A  (+  b  sin  0  +  -  y)] 

u  x  y 


Using  the  expressions  for  A  ,  A  ,  0+  in  terms  of 

x  y 


simplifies  this  to 


_  _g_ 


±  4U2  ®+ 


To  lowest  order  (in  0  )  we  have 


.2  32A 

3  4  ~  gx  x 

2*2  2  ' 

3A  U  3A 

y  y 


32A 


But 


-2  0. 


3  A 


-i . .  „  0.2  34 

“>  T(A  -A  )  - = 

,  7  v  v  7  7 


—  4  TT  /  A 


The  results  for  1+  are  then 


I 


+ 


2  —  - 
2  iTU  +14  -14 

-  e~  e 

g 


/ 


b(-x 


x  e  g  (0.)  h  (0+) 

4U  0+  - 


To  recapitulate: 


0 


2  0_a 


(Remember  we  have  assumed  y  >  b  ! ) 


« \r> 


V.  Explicit  Forms 


A.  The  Contributions  From  The  Stationary  Point  Where 

l  ~  l  -1,1  -  0 

x  y 

Here  0  =  -  y/x  . 

and  0+  *  0  ,  0  *  it  -  0 


From  the  results  of  Section  IV,  it  is  readily  seen  that  all  of  1^  +  I2 

2 

+  vanish  at  least  as  0  . 


is 


The  only  term  in  I^  +  1^  +  1^  which  does  not  vanish  similarly 


d*  i 

y  x 

£(1-2A) 


-i  -f  x  +  l  y) 

u2  x  y 


7T 

x  /  cos  9  d  9 
o 


[e+  +  e_] 


+  c.c. 


Using  equation  (46)  we  then  find 


-i*  (x  - 

f  4  IT 
i  e  e 


-1-f  [<*«*>- j’l 


x  '  -  +ii  V 


B.  Contributions  From  the  Far  Stationary  Point 


Here 


40* 


40. 


20. 


•+  ■  7  -  >  ^ 


1+2  S-t 

2  £  b 


0  »  . 

+  x 


Then 


sin  0+ 


1  +  0  (02) 


GH-a 


cos  9+  =*  +  2 


+  0  (0p 


One  might  suspect  that  because  of  the  overall  factor  of 
(0+)  ^2  in  equation  (47)  that  there  may  be  singularities  in  some 
derivatives  of  our  potential.  However,  it  can  be  shown  that  almost 
all  the  terms  in  equations  (30)  and  (31)  are  at  worst  of  order 


We  consider  two  examples: 


20+a  (  20+  ¥ 

b  1  20. 


-  ^  |  e+  »  0  to  0  (03) 


As  a  consequence  we  see  that  the  dominant  contrilbutlon  for 


0  <  <  1  is: 


dl  i  ~i  [V  +  ly  yl 

iUb  »  y  x  U 
2u  J-»  1-2*  6 


x  J  co s  6  d6  5  f  (£  )  t-V1  ] —  [e+  +  e  ] 
‘  o  o  ,.2  2.  1 

o  (G  ~W  ) 


+  c.c. 


Using  our  rules  for  integration  this  becomes 


<j>  »  -2Ua 


f«o> 


g  b(-x)  c 


3  ,  7 

r  -2  +102 

X  *  °+  6  4U20+ 


-  ig* 
4U2  0 


+  c.c. 


VI. 


Discussion 


We  conclude  that  with  our  approximations  the  dominant  term 
for  the  potential  near  the  axis  is  given  by  equation  (52).  This  has  a 
singularity  at  y  »  b  not  at  y  ■  0. 

The  potential  and  displacement  at  the  singularity  are  both 
zero.  But  the  slopes  are  infinite. 

Where  does  the  singularity  come  from?  We  think  not  from  our 
2 

use  of  the  approximation  U  /gL  <  <  1  .  Indeed  tracing  back  to  the 
origin  of  the  term  we  see  it  arises  from  3$/3x  on  the  waterline.  Our 
procedure  guarantees  that  3$/3x  there  la  represented  exactly. 

Rather  the  origin  would  seem  to  result  from  an  inappropriate 
use  of  the  stationary  phase  method  when  evaluating  equation  (51). 

(This  does  not  seem  to  matter  for  all  other  terms  of  equations  (30) 
and  (31).  They  are  multiplied  by  a  high  powers  of  the  vanishing 
quantity.)  The  problem  with  evaluating  equation  (51)  by  stationary 

phase  is  that  - j  *  0  as  0  ♦  0 

3 1  2 

y 

Normally  one  would  handle  such  a  situation  by  proceeding  to 
the  next  term  in  the  Taylor  series  of  ♦  around  the  stationary 


point.  Thus  instead  of  the  Fresnel  type  integral  one  would  have  an 
Airy  integral  to  describe  the  result  in  the  vicinity  of  the  bogus 

singular  point.  Here,  however,  one  readily  shows  that  not  only  does 

2  3 

j*1  *  3  ^  4 

- r-  +  0  but  so  does  - ^  .  Indeed  all  higher  derivatives  go  to  zero 

3 1  1  3 1 

as  y0+  +  0  .  (They  evenygo  to  zero  faster  and  faster  as  the 
derivative  increases.) 

We  conclude: 

(1)  The  result  of  equation  (52)  would  seem  to  indicate  that 
there  are  large  slopes  in  the  wake  near  but  not  on  the  axis. 

(2)  A  definitive  conclusion  (of  the  linear  theory)  awaits  a 
better  evaluation  of  the  integrals  of  equation  (51). 

(3)  Since  we  have  used  a  very  specific  hull  model  it  is  not 
immediately  obvious  as  to  what  the  shape  dependence  really  is. 

However,  the  procedure  outlined  here  strongly  suggests  that  large 
slopes  should  be  produced  in  the  wake  at  distances  from  the  axis 
comparable  to  some  measure  of  the  ship  width. 
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